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Abstract 



We report a new measurement of the pseudorapidity (r/) and transverse- 
energy (Et) dependence of the inchisivc jet production cross section in pp 
colhsions at ^/s = 1.8 TeV using 95 pb~^ of data collected with the D0 
detector at the Fermilab Tevatron. The differential cross section d'^a/ {dExdri) 
is presented up to |r7| = 3, significantly extending previous measurements. The 
results arc in good overall agreement with next-to-leading order predictions 
from QCD and indicate a preference for certain parton distribution functions. 
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This past decade has witnessed impressive progress in both the theoretical and exper- 
imental understanding of the coUimated streams of particles or "jets" that emerge from 
inelastic hadron collisions. Theoretically, jet production in hadron collisions is understood 
within the framework of Quantum Chromodynamics (QCD), as a hard scattering of the con- 
stituent partons (quarks and gluons) that, having undergone a collision, manifest themselves 
as jets in the final state. QCD predicts the amplitudes for the hard scattering of partons 
at high energies. Perturbative QCD calculations of jet cross sections using accurately 

determined parton distribution functions (PDFs) @J^, have increased the interest in jet 
measurements at the ^/s = 1.8 TeV Tevatron proton-antiproton collider. Consequently, the 
two Tevatron experiments, D0 and CDF, have served as prominent arenas for studying 
hadronic jets. 

In this Letter, we report a new measurement of the pseudorapidity (77) and transverse- 
energy {Et) dependence of the inclusive jet production cross section which examines the 
short-range behavior of QCD, the structure of the proton in terms of PDFs, and possible 
substructure of quarks and gluons. We present the differential cross section d'^a/ {dETclrj) as 
a function of jet Et in five intervals of r], up to \ri\ = 3, where the pseudorapidity is defined 
as 7] = In [cot (^/2)], with 6 being the polar angle. The present measurement is based 
on 95 pb~^ of data collected with the D0 detector |^ during 1994-1995, and significantly 
extends previous measurements P], as indicated by the kinematic reach shown in Fig. |l[ 

The primary tool used for jet detection is the compensating, finely segmented, liquid- 
argon/uranium calorimeter, which provides nearly full solid-angle coverage {\ri\ < 4.1). Jets 
are defined and reconstructed off-line using an iterative fixed-cone algorithm with a cone 
radius ofTl = 0.7 in the rj-ip space, where (p is the azimuth. The missing transverse energy 
{^t) is calculated from a vector sum of the individual Et values in all the cells of the 
calorimeter. Calorimeter cells can occasionally provide spurious noise signals; to diminish 
their effect on jets, such cells are identified and suppressed using specific on-line and off-line 
algorithms. 

During data taking, events were selected with a multi-stage trigger system. The first 
stage signaled an inelastic pp collision. In the next stage, the trigger required a jet in a 
calorimeter region of At] x A(f = 0.8 x 1.6, with Et above a preset threshold. In the 
last trigger stage, selected events were digitized and sent to an array of processors. Jet 
candidates were reconstructed using a cone algorithm, and the entire event was recorded if 
any jet Et exceeded a specified threshold. The four software filters used in this analysis had 
Et thresholds of 30, 50, 85, and 115 GeV, and accumulated integrated luminosities of 0.364, 
4.84, 56.5, and 94.9 pb"^ respectively p. To present the full range of the data, the cross 
sections obtained from the four jet filters are combined in contiguous regions of Et in such 
a way that the more restrictive trigger is adopted as soon as it is more than 99% efficient. 

The position of the primary interaction vertex is reconstructed using data from the central 
tracking system. The two vertices with the largest number of associated tracks are retained 
for further analysis. At high instantaneous luminosities, multiple interactions are common, 
and to correct for inefficiency of the tracking system in identifying the primary vertex, we 
use the global event quantity 5't = | Et^ \ . The vertex with the smaller value of St is 
defined as the correct event vertex, and all kinematic variables are calculated with respect 
to it. The dependence of jet Et on luminosity was studied, and found to be negligible. 
At high pseudorapidities, the jet reconstruction algorithm introduces a bias towards 77 = 
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FIG. 1. The kinematic reach of this measurement along with that of other cohider and 
fixed-target experiments in the plane of the parton momentum fraction x and the square of the 
momentum transfer Q^. 



0. Furthermore, the Snowmass jet reconstruction algorithm used in the theoretical 
predictions has a different definition for jet angles than that used in the standard D0 
off-line algorithm. Jet rj values are corrected for this difference, which also removes any 
instrumental bias in reconstruction of jet polar angles 0. 

Backgrounds introduced by electrons, photons, detector noise, accelerator losses, or cos- 
mic rays are removed using quality criteria developed for jets with Ir^l < 3. To preserve 
the pseudo-project ive nature of the D0 calorimeters, the longitudinal {z) position of the 
interaction vertex is required to be within 50 cm of the detector center; this requirement is 
(88.7 ± 0.1)% efficient. A cutoff on 

removes background from cosmic ray showers and 
misvertexed events. 4t must be smaller than the lesser of 30 GeV or O.SEjn of the leading jet 
if the leading jet is central (|?7| < 0.7), or less than 0.7Et otherwise. This criterion is nearly 
100% efficient. Jet quality is based on the pattern of energy deposition in the calorimeter. 
The combined efficiency for jet quality ranges from about 99.5% at lowest Erp and |?7| to 
approximately 98% at highest Ft and |?7|. 

The jet energy calibration, applied on a jet by jet basis, corrects (on average) the re- 
constructed Et for variation in the hadronic response of the calorimeter, for the energy 
associated with underlying spectator interactions, for multiple pp interactions in the same 
crossing, noise originating from uranium decay, the fraction of any particle's energy that 
showers outside of the reconstruction cone, and for detector non-uniformities. A complete 
discussion of the jet energy calibration can be found in Ref. |Tl]]. An independent test of the 
jet energy scale, based on the balance in transverse energy in photon-jet and jet-jet data, 
confirms the validity of the D0 jet-calibration procedure up to |?7| = 3 ||^. 

In each bin of tj-Et, the average differential cross section, d'^a/{dETdr]), is calculated as 
N/{AriAET eJCdt), where AtjAEt is the tj-Et bin size, N is the number of jets observed 
in a bin, e is the total overall efficiency for jet and event selection, and JCdt represents the 
integrated luminosity of the data sample. Statistical uncertainties in the values of the cross 
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sections are defined by one standard deviation Poisson fluctuations in the associated N. 

Energy resolution of tlie D0 calorimeters distorts the jet cross section in Ex. Although 
the resolution is essentially Gaussian, the jet cross section is shifted to larger Et due to the 
steeply falling dependence of jet production on Et- This effect is removed from the data 
through an unfolding procedure. We measure the fractional jet energy resolutions based on 
the "same side" {rirri2 > 0) subset of dijet events in the data sample. Using the imbalance 
in the Et of the two leading jets, in each interval of \r)\, we parameterize the fractional 
jet energy resolution as a function of jet Ej-, following the standard description of single- 
particle energy resolution, based on the noise, sampling, and constant terms. To determine 
the amount of distortion in the cross section in each of the five |?7| intervals, we take an 
ansatz function of the form e'^E^{l + '~)2Et/ ^/s) , numerically smear it according to the 
parameterized resolution in each E^ bin, and fit this smeared hypothesis to the observed 
cross sections to extract the five sets of four free parameters, a, /3, 7, and 5. The bin-by- 
bin ratio of the original over the smeared ansatz for each range of |?7| gives the unfolding 
correction with which we rescale the observed cross section to remove the distortion from 
jet energy resolution 0. 

The jet angular resolution is very good at all 77, and its effect on the cross section is 
negligible, but it is possible to distort the jet polar angle through a mismeasurement of the 
z position of the vertex. However, a Monte Carlo study demonstrates that such effects are 
negligible because distortions in jet Et are nearly fully compensated by bin-to-bin migrations 
in |?7| from the smearing of the z coordinate of the vertex 0. 

The final measurements in each of the five \ri\ regions, along with statistical uncertainties, 
are presented in Fig. || (tables of the measured cross sections can be found in Refs. ||,|2 



The measurement spans about seven orders of magnitude and extends to the highest jet 
energies ever reached. Figure ^also shows 0{al) theoretical predictions from JETRAD 
with renormalization and factorization scales set to half of the Et of the leading jet and 
using the CTEQ4M PDF. 

Figures |] and § provide more detailed comparisons to predictions on a linear scale for 
several PDFs (for other PDFs, see Ref. 0). The error bars are statistical, while the shaded 
bands indicate one standard deviation systematic uncertainties. Because the theoretical 
uncertainties due to variations in input parameters are comparable to the systematic uncer- 
tainties |T^, these qualitative comparisons indicate that the predictions are in reasonable 



agreement with the data for all \rj\ intervals. 

To quantify the comparisons, we employ a specially derived and previously studied 
statistic of the form = T.i,j iPi - T^) [(T^/Di) C^,- (T./D,)]'^ (D,- - T,), where 

(Dj — Tj) is the deviation of the measured cross section (Dj) from the prediction (Tj) in the 
z-th bin, Cjj is the full covariance matrix of the measurement defined as YliaP%^'i^'j ■> 
where the sum runs over all sources of uncertainties, pij is the correlation coefficient be- 
tween the i-th and j-th bins, and ai is the uncertainty in the i-th bin. The T/D factors are 
introduced to reduce the bias towards lower values of originating from highly correlated 
systematic uncertainties present in Cy 0. There are 90 t]-Et bins in this measurement. 

While the statistical uncertainties are not correlated in Et or 77, the systematic uncer- 
tainties are fully correlated in both variables except for: (i) efficiencies for data selection, 
which are uncorrelated in 77, (ii) parameterizations of jet energy resolutions and fits to the 
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FIG. 2. The single inclusive jet production cross section as a function of jet Et, in five pseu- 
dorapidity intervals, showing only statistical uncertainties, along with theoretical predictions. 
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FIG. 3. Comparisons between the D0 single inclusive jet cross sections and the O(a^) QCD 
predictions calculated by JETRAD with the CTEQ4HJ (•) and CTEQ4M (o) PDFs. The highest 
Et points are offset slightly for CTEQ4M. 



unfolding ansatz, which are uncorrelated in t], (ill) the hadronic response, which is partially 
correlated in Et and rj, with the correlation matrix in terms of average bin energies given 
in Ref . []TT| . Uncertainties in the showering correction arise dominantly from the lack of full 
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FIG. 4. Comparisons between the D0 single inclusive jet cross sections and the 0(a'^) QCD 
predictions calculated by JETRAD with the MRSTgt (•) and MRST (o) PDFs. The highest Et 
points are offset slightly for MRST. 



agreement of the lateral shower profiles observed in the data and in the Monte Carlo. The 
residual discrepancy is similar for all Et and rj regions. Consequently, the correlations of 
the showering correction are large in Et 0] as well as in rj. Uncertainties due to jet en- 
ergy calibration are the dominant source of error in the cross section and range from about 
12-20% at lowest Et to about 35-80% at highest Et, getting larger with rj for a fixed Et- 
They are driven by the uncertainties due to the hadronic response parameterization at high 
Et and due to the showering correction at high Et and, notably, at high rj. The second 
largest source of uncertainty is the jet energy resolution parameterization and the unfolding 
procedure which typically gets worse at low and at high Et and ranges from about 3-5% at 
lowest Et to about 10-20% at highest Et- These are followed by the uncertainties due to 
integrated luminosity which are approximately 6% (8%) for the data collected with the jet 
filters with two highest (lowest) Et thresholds, and by the uncertainties due to data selection 
which are on the order of 1% throughout the dynamic range of the measurement 0. 

For all PDFs we have considered. Table | lists the x^y X^/dof, and the corresponding 
probabilities for 90 degrees of freedom (dof). We have verified that the variations of cor- 
relation coefficients within the range of their uncertainties give a similar ordering of the 
X^, hence a similar relative preference of PDFs. The absolute values of and associated 
probabilities vary somewhat with variations in the correlations in Et and, to a much lesser 
extent, with variations of correlations in t]. The theoretical predictions are in good quan- 
titative agreement with the experimental results. The data indicate a preference for the 
CTEQ4HJ, MRSTgt, and CTEQ4M PDFs. The CTEQ4HJ PDF has enhanced gluon con- 
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tent at large x, favored by previous measurements of inclusive jet cross sections at small 



r] 115|,|1|], relative to the CTEQ4M PDF. The MRSTgt PDF includes no intrinsic parton 
transverse momentum and therefore has effectively increased gluon distributions at all x 
relative to the MRST PDF. 

TABLE I. The x^i X^/dof, and the corresponding probabihties for 90 degrees of freedom for 
various PDFs studied. 



PDF 




xVdof 


Probability 


CTEQ3M 


121.56 


1.35 


0.01 


CTEQ4M 


92.46 


1.03 


0.41 


CTEQ4HJ 


59.38 


0.66 


0.99 


MRST 


113.78 


1.26 


0.05 


MRSTgi 


155.52 


1.73 


<0.01 


MRSTgt 


85.09 


0.95 


0.63 



In conclusion, we have reported a new measurement of the pseudorapidity and transverse- 
energy dependence of the inclusive jet cross section in proton-antiproton collisions at 
^/s = 1.8 TeV. Our results extend significantly the kinematic reach of previous studies, are 
consistent with QCD calculations over the large dynamic range accessible to D0 {\ri\ < 3), 
and indicate a preference for certain PDFs. Once incorporated into revised modern PDFs, 
these measurements will greatly improve our understanding of the structure of the proton 
at large x and Q^. 
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